Gamma-rays from dark matter annihilations strongly constrain the substructure in 

halos 
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Recently, it has been shown that electrons and positrons from dark matter (DM) annihilations 
provide an excellent ht to the Fermi, PAMELA, and HESS data. Using this DM model, which 
requires an enhancement of the annihilation cross section over its standard value to match the 
observations, we show that it immediately implies an observable level of 7-ray emission for the 
Fermi telescope from nearby galaxy clusters such as Virgo and Fornax. We show that this DM 
model implies a peculiar feature from final state radiation that is a distinctive signature of DM. 
Using the EGRET upper limit on the 7-ray emission from Virgo, we constrain the minimum mass 
of substructures within DM halos to be > 5 x 10~ 3 Mq - about four orders of magnitudes larger 
than the expectation for cold dark matter. This limits the cutoff scale in the linear matter power 
spectrum to k < 35 kpc -1 which can be explained by e.g., warm dark matter. Very near future Fermi 
observations will strongly constrain the minimum mass to be > 10 3 Mq : if the true substructure 
cutoff is much smaller than this, the DM interpretation of the Fermi/PAMELA/HESS data must 
be wrong. To address the problem of astrophysical foregrounds, we performed high-resolution, 
cosmological simulations of galaxy clusters that include realistic cosmic ray (CR) physics. We 
compute the dominating 7-ray emission signal resulting from hadronic CR interactions and find that 
it follows a universal spectrum and spatial distribution. If we neglect the anomalous enhancement 
factor and assume standard values for the cross section and minimum subhalo mass, the same model 
of DM predicts comparable levels of the 7-ray emission from DM annihilations and CR interactions. 
This suggests that spectral subtraction techniques could be applied to detect the annihilation signal. 

PACS numbers: 95.35.+d, 95.85.Pw, 98.62.Gq, 98.65.-r, 98.70.Sa 



The identity of the DM in the Universe has been the 
subject of intense speculation. In particular, the hierar- 
chical formation of structure, as indicated from numer- 
ical simulations of Cold DM (CDM), agrees very well 
with observations on scales of galaxy clusters and larger, 
whereas the small scale behaviour on galactic and sub- 
galactic scales is more unsecure. If dark matter particles 
have weak interactions, one would expect possible signals 
from annihilations (or decays). 

Data from a new generation of cosmic ray detectors 
have indeed been tentatively interpreted in terms of such 
signatures of DM. In particular, the positron fraction 
measured by the PAMELA satellite [lj and the sum of 
electrons and positrons by ATIC Q have shown an un- 
expected excess. Very recent data from Fermi-LAT 0, Q 
and H.E.S.S. [f| on the sum of electrons and positrons do 
not confirm the peak claimed by ATIC, but still indicate 
an excess compared to the expected background in con- 
ventional models. A number of attempts have been made 
trying to explain the excess due to DM annihilation (for 
recent reviews, see d, 0) but also other astrophysical 
sources such as pulsars have been investigated (e.g., 0, 
and references therein). 

One class of DM models that fits the new data has halo 
annihilation primarily into muonpairs, which then decay 
to electrons and positrons. In |9[ examples of fits with 



H.E.S.S. data) were obtained by such DM models. It was 
pointed out that if the annihilation goes directly into a 
pair, a striking signature may be present. This 
is caused by the direct emission of 7-rays from the hnal 
state (final state radiation, FSR), which gives a peculiar 
energy spectrum, with E 2 dN 1 /dE almost linearly rising 
with energy The same, but weaker, feature may in fact 
exist also for the theoretically perhaps more easily moti- 
vated models with intermediate spin-0 boson decay, but 
in this Letter we only treat the somewhat simpler direct 
annihilation case. In this Letter we show that such DM 
models face considerable tension with existing 7-ray lim- 
its from clusters of galaxies, systems which will be very 
interesting to detect and study with coming 7-ray detec- 
tors (Fermi, H.E.S.S., MAGIC, VERITAS, and eventu- 
ally CTA). 

Galaxy clusters constitute the most massive objects in 
our Universe that are forming today. This causes their 
DM subhalo mass function to be less affected by tidal 
stripping compared to galaxy sized halos that formed 
long ago. The annihilation luminosity of the smooth DM 
halo component scales as 
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remarkable quality (which also fit PAMELA [1| and new where the virial mass M200 and the concentration c (see 
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Eqn. [3]) arc the two characteristic parameters of the uni- 
versal Navarro- Frenk- White (NFW) density profile pnfw 
of DM halos (lQj . Hence, the flux ratio of a nearby clus- 
ter (Virgo) to a prominent dwarf spheroidal (Draco) is 
given by 



F. 



cluster 



Fdv 



80kpc \ 2 /2 x 1O 14 M 
17MpcJ ^ 10 s M, 



i 



(2) 



assuming an early formation epoch of the dwarf galaxy 
before the end of reionization [11| . Once a satellite galaxy 
is accreted by our Galaxy, the outer regions are severely 
affected by tidal stripping. The longer a satellite has 
been part of our Galaxy, and the closer it comes to the 
center during its pericentral passage, the more material 
is removed [12j. In contrast, the substructure in clus- 
ters is not affected in the outer regions and may en- 
hance the DM annihilation signal over its smooth con- 
tribution considerably as we will see in the following. 
The FSR feature of DM annihilation may in addition 
be more easily visible in clusters, as the average inten- 
sity of starlight, which may give a masking signal due 
to inverse Compton scattering of the copiously produced 
electrons and positrons, is lower than in the Milky Way 
MW. For previous work related to dark matter in clus- 
ters, see, e.g., 

[HEi]. All halo masses and length scales 



are scaled to the currently favored value of Hubble's con- 
stant, Ho = 70 km s _1 Mpc - . We define the virial mass 
M200 and virial radius r2oo as the mass and radius of a 
sphere enclosing a mean density that is 200 times the 
critical density of the Universe p CI . 

Method. As our default model for dark matter an- 
nihilation, we take the Sommerfeld-enhanced (see, e.g., 
[TBI , [lij ) direct muon annihilation mode of Q, i.e., 
mass m x = 1600 GeV and effective enhancement factor 
1100 relative to the standard annihilation cross section 
(av)o ~ 3x 10 -26 cm 3 s" 1 . It is non-trivial to rescale this 
boost to the corresponding value for the cluster environ- 
ment. This may either give a smaller or larger value, 
depending, e.g., on the velocity dispersion of bound sub- 
structure in the cluster, and whether the Sommerfeld en- 
hancement (SFE) increases down to very small velocities, 
or if it saturates at some minimum velocity [l7j . We 
choose a simple and generic model for the SFE factor 
Bsfe(i7«) « 0.7 c/a v and saturate at 

&v. min — ■ 

200 km/s 

[3^ . For a given cluster, we assume a constant vel ocity 
dispersion of a v = 960 km/s x (M 200 /10 15 Mq) 1 / 3 [Hf. 
This scaling would result in a boost factor of B s f c = 220 
for a cluster with M200 = 10 15 Mq. 

Using cluster masses from the complete sample of the 
X-ray brightest clusters (the extended HIFLUGCS cat- 
alogue, [19j) we identify the brightest clusters for DM 
annihilation. In models with SFE, the DM flux to lead- 
ing order scales as a power- law F ~ M 20 ]/ 3 M^ ^ /D 2 . 
The first factor accounts for the SFE and the second one 
is derived from Eqn. [TJ using a power-law fit to the mass 



dependence of the NFW halo concentration derived from 
cosmological simulations with M200 ^ 1O 1O M0 [20l |. 
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Note that Eqn.[3]agrees well with [2l| for cluster-mass ha- 
los after converting the concentration definitions accord- 
ing to jH). This yields Fornax[M 2 oo = 10 14 M Q ) and 
Virgo (M200 = 2.1 x 10 14 M o , ^) as the prime targets 
for DM observations and we additionally decide in favor 
of the well studied cluster Coma (M 20 o = 1.4 x 10 15 M o ) 
for comparison. 

The differential photon flux from DM annihilation 
within a given solid angle Af2 along a line-of-sight (los) 
is given by 
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where <? S m {E-y, r) is the source function from the smooth 
halo with contributions from two main processes: DM 
annihilating to p + / p~ which decay to e + /e~ pairs that 
Compton up-scatter CMB photons (IC) and FSR. The 
source function of FSR is given by 
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where the annihilation rate density = 
(p/m x ) 2 (av)i/2. The i runs over all 7-ray pro- 
ducing channels each with the spectrum d j^7'' and 
annihilation cross-section (eru);. We use the standard 
photon distribution for the final state radiation from 
our DM model annihilating directly to charged leptons 
assuming m x 3> mi [24| . For the remaining part of this 
work, the Einasto density profile for DM halos j2f| is 
used, normalized with po = Pnfw(?"s)/4 relying on the 
assumption that 90% of the flux from a NFW density 
profile and an Einasto density profile originate from 
within the scale radius r s = r2oo / c 

The product of enhancement factors from SFE B s f c (a v ) 
and from substructure enhancement over the smooth halo 
contribution B su b (r) = 1 + q su b( r ) I Qsm( r ) is denoted by 
Bp(a v ,r) = B s f e ((T t ,) B su b(r). High-resolution DM sim- 
ulations of the MW suggest an enhancement from sub- 
structures of approximately 220 inside r2oo assuming that 
the subhalos extrapolate smoothly down from the simu- 
lation resolution limit to smallest scales [2(|, with most 
of the substructure residing in the outer part of the MW 
halo. We fit the luminosity L su b = / d-E^j dVq su t, from 
substructures inside a radius r following [271 ] . 

L sub (< r) = 0.8CL sm (r 200 ) (r/r 2 oo) - 8(r/r20or °' 315 , 

(6) 

where L sm (r2oo) is the smooth cluster halo luminosity 
within r2oo- The normalization C = (M m i n /Afii m ) ' 226 , 
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where Af m i n = 10 Mq is the minimum subhalo mass 
in the simulation and M\\ m the free streaming mass. 
While its conventional value is 10 -6 Mq (28|, we will 
constrain this quantity by requiring consistency with the 
non-detection of 7-ray emission from clusters by EGRET: 
the smaller Mn m , the more substructure is present and 
the larger is the expected 7-ray signal. This approach of 
fitting the scaling behaviour of L su b(A/n m ) directly from 
numerical simulations self-consistently accounts for the 
radial dependence of the substructure concentration [26| . 
We note that this might result in a slight overestimate 
of the substructure luminosity if the assumed power-law 
scaling flattens towards smaller scales although current 
simulations show no sign of such a behavior which is also 
not expected since we are approaching the asymptotic 
behavior in the power spectrum on these scales. 

The source function of inverse Compton emission re- 
sulting from DM annihilating is given by 
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where Pic is derived by convolving the IC cross-section 
with the differential target photon number density 
Assuming that the spatial diffusion time scale is much 
larger than the energy loss time scale, the total equilib- 
rium distribution of the electrons plus positrons is given 
by 
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where ^ L denote the differential number of electrons 
plus positrons resulting from an annihilation event, 
-Bcmb = 3.24(1 + z) 2 fiG denotes the equivalent field 
strength of the CMB, and we parametrize the magnetic 
field of the galaxy cluster by B(r) = 3fiG [n e (r)/n e (0)]°- 7 , 
which follows from flux frozen magnetic fields. 

To address the problem of source confusion by astro- 
physical backgrounds, we perform high-resolution, cosmo- 
logical simulations [29| of a sample of 14 galaxy clusters 
[30J . They span over one and a half decades in mass and 
follow radiative gas physics, star formation, supernova 
feedback. In particular, we use an updated version of 
the cosmic ray physics that is capable of following the 
spectral evolution of the cosmic ray (CR) distribution 
function by tracking multiple CR populations - each be- 
ing described by its characteristic power-law distribution 
with a distinctive slope that is determined by the accel- 
eration process [3l|. We compute the dominating 7-ray 
emission signal from decaying neutral pions that result 
from hadronic CR interactions with the ambient gas fol- 
lowing We find that it obeys a universal spectrum 
and spatial distribution (details will be presented in [3l|). 
This allows us to reliably model the cosmic ray signal 



from nearby galaxy clusters using their true density pro- 
files as obtained by X-ray measurements [33| that we map 
onto our simulated density profiles. We compute 7-ray 
luminosity-mass scaling relations of our sample [34| and 
use these to normalize the CR induced emission of all 
clusters in HIFLUGCS 0. In our optimistic CR model, 
we calculate the cluster's total 7-ray flux within a given 
solid angle while we cut the emission from our individual 
galaxies and compact galactic-sized objects in our more 
conservative baseline model 



Results and Discussion. In Fig. [TJ we find that, 
given our assumptions, the DM annihilation signal in 
Fornax and Virgo should be clearly visible by Fermi. It 
dominates over the CR induced signal for both of our CR 
models. The annihilation signal within 3.5 deg is similar 
for Fornax and Virgo - these are clusters with a compa- 
rable distance but the latter being twice as heavy: the 
larger signal of Virgo due to its larger mass is counter- 
acted by the larger substructure boost of Fornax for the 
same angular extent. Using the standard assumptions 
for the limiting mass of substructures within DM halos 
of 10~ 6 M Q , we show that the resulting annihilation flux 
from Virgo is already in conflict with the EGRET upper 
limit. This allows us to place a lower bound on the limit- 
ing mass Mum = 5 x 10 -3 M Q and hence to constrain the 
free streaming scale in the linear matter power spectrum 
to 
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The Fermi sensitivity will allow us to place an even more 

El, (9) 

stringent limit of Mu m > 1O 3 M which is approaching 



the upper limit M\ im < 2 x 10 M Q derived from Ly-a 
power spectrum measurements 35|. The contribution of 
the smooth DM halo component towards high galactic 
latitudes within 3.5 degree amounts to F{> 100 MeV) ~ 
9 x 10 _9 ph cm _2 s _1 which should be easily detectable 
by Fermi, especially considering an enhancement factor 
of a few from substructure in the MW halo. This al- 
lows us to finally scrutinize the DM models motivated by 
the recent Fermi, H.E.S.S., and PAMELA data. Assum- 
ing SFE, the next generation of imaging air Cerenkov 
telescopes have good prospects of detecting the DM an- 
nihilation signal while it will be very difficult without an 
enhancement. We show in Fig. [21 that the DM annihila- 
tion flux is substantially boosted due to substructures in 
clusters as well as in the MWs halo that has a smooth an- 
gular emission characteristic but is expected to have the 
same spectral behavior. This provides hope that even in 
the absence of SFE, the DM annihilation flux is of the 
same order of magnitude as our conservative model of 
CR induced 7-ray emission. The very distinctive spectral 
properties of the DM-induced 7-rays and the universality 
of the CR spectra suggest that spectral subtraction tech- 
niques could be applied to detect the annihilation signal 
and characterize the properties of DM. 



4 




FIG. 1: Left: differential spectra for 3 different clusters observed within a solid angle of diameter A9 — 3.5 degree 
(neglecting the contribution of the MW's halo). We show the inverse Compton up-scattered CMB photons from e + /e~ 
pairs that themselves result from DM annihilating to fj, + 1 pT (solid) as well as the final state radiation (dash-dotted). We 
include Sommerfeld enhancement (SFE) and the enhancement from cluster substructures down to a limiting substructure 
mass of Mum = 5 x ICP 3 M©. The 7-ray emission from decaying neutral pions that result from hadronic cosmic ray (CR) 
interactions with the ambient gas is shown with dashed lines for our conservative model without galaxies. Right: we 
contrast the integrated spectrum of Virgo for the EGRET angular resolution, AO = 5.8 (red), with that of imaging air 
Cerenkov telescopes, AO — 0.1 (blue), and compare those to the point source sensitivity curves on the 5<r level of Fermi (2 
year all-sky survey) as well as MAGIC2 and CTA (50 hours). We choose M\i m = 5x 10~ 3 M©, so that the resulting flux 
is just consistent with the EGRET upper limit [3(|. This reduces the substructure boost from ~ 220 to 50. 




FIG. 2: Studying separately the contribution from substructure and the SFE to the differential spectra of Fornax within 
a solid angle of A8 = 3.5 degree. Left: the dark matter annihilation spectra with SFE and Mn m = 5 x 10 -3 Mq (red) 
is compared to the pion decay spectrum from CR interactions (blue). We show the pure contribution from the smooth 
cluster halo (dotted, B s f c = 490,B su b = 1), to which we add substructures (dashed, B s f c = 490, B su b = 50), and to which 
we additionally add the line-of sight contribution due to the Milky Way's halo towards Fornax (solid). Right: we compare 
the hadronically induced pion decay spectrum to the DM annihilation signal without SFE. The pion decay spectrum shown 
with (dotted) and without (solid) the contribution due to galaxies and dense point sources. The substructure boosted 
DM annihilation signal including the MW contribution (dash-dotted), but assuming a standard value for the limiting 
substructure mass of Mn m = 10 _6 M©. 



Conclusions. The DM models motivated by the re- 
cent Fermi, H.E.S.S., and PAMELA measurements re- 
quire an anomalous boost factor of 1100. Assuming that 
SFE entierly accounts for this boost, this necessarily pre- 



dicts large annihilation fluxes from nearby galaxy clus- 
ters even in the case of somewhat reduced SFE due to 
the larger velocity dispersion of clusters. Using standard 
assumptions for the limiting mass of substructures within 
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DM halos, we find a violation of the EGRET upper limit 
in Virgo. The lighter a DM particle, the larger the in- 
duced free-streaming scale in the power spectrum and 
the higher the mass cutoff for the smallest structures: 
since the DM interpretation fixes the DM particle mass 
this locks in a minimum substructure mass. Hence, a 
non-detection of 7-rays at the predicted level by Fermi 
would provide a serious challenge for the standard as- 
sumptions of the CDM power spectrum, or call for a new 
dynamical effect during non-linear structure formation 
that wipes out the smallest structures. The resolution 
may of course also be that the rising positron ratio mea- 
sured by PAMELA and the electron plus positron excess 
seen by H.E.S.S. and Fermi is caused by local astrophysi- 
cal sources, e.g. pulsars, and is unrelated to dark matter. 
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